Multiphoton transitions in a macroscopic quantum two-state system 
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We have observed multiphoton transitions between two macroscopic quantum-mechanical super- 
position states formed by two opposite circulating currents in a superconducting loop with three 
Josephson junctions. Resonant peaks and dips of up to three-photon transitions were observed in 
spectroscopic measurements when the system was irradiated with a strong RF-photon field. The 
widths of the multiphoton absorption dips are shown to scale with the Bessel functions in agreement 
with theoretical predictions derived from the Bloch equation or from a spin-boson model. 
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A macroscopic quantum two-state system (TSS) of- 
fers a unique testing ground for exploring the founda- 
tions of quantum mechanics Q . This system can be in 
a quantum-mechanical superposition of macroscopically 
distinct states, and its quantum nature can be revealed 
by measuring the absorption of an integer number of pho- 
tons from an externally applied photon field 2]. Since a 
macroscopic quantum system cannot be completely de- 
coupled from its environment, dissipative and decoher- 
ence effects are unavoidable jli, ■ In addition to the 
investigation of fundamental physics, the quantum TSS 
also serves as an elementary carrier of information in a 
quantum information processor in the form of a quantum 
bit (qubit) Q . Artificially designed superconducting cir- 
cuits with mesoscopic Josephson junctions constitute an 
important class of macroscopic quantum systems. The 
charge degree of freedom of Cooper pairs is used to induce 
coherent quantum oscillations between two charge states 
of a Cooper pair box Q . A circuit with a single relatively 
large Josephson junction, which is current-biased close 
to its critical current, forms a so-called Josephson phase 
qubit 6]. Moreover, three Josephson junctions arranged 
in a superconducting loop threaded by an externally ap- 
plied magnetic flux constitute a flux qubit The de- 
vice could be prepared in a quantum superposition of two 
states carrying opposite macroscopic persistent currents 
jg^. Coherent Rabi oscillations have been reported, when 
the qubit and the readout device are connected to ob- 
tain a large signal 0. Rabi oscillations have also been 
observed in a system, where the qubit and the readout 
device are spatially separated |0| . Since these solid-state 
devices are thought to be scalable up to a large number 
of qubits, they are of particular interest in the context of 
solid-state quantum information processing 11]. 

The energy scale of quantum circuits containing 
Josephson junctions is in the microwave regime. This 
property was demonstrated in the current-voltage char- 
acteristics of a Josephson junction under microwave ir- 
radiation displaying the well-known Shapiro steps [l^ . 
They appear at voltages corresponding to integer multi- 



ples of the applied microwave energy. With this phe- 
nomenon, the superconductor phase difference at the 
junction, which is a macroscopic degree of freedom, can 
be treated as a classical degree of freedom moving in 
the Josephson potential. In contrast, the quantum- 
mechanical behavior of the macroscopic phase difference 
was demonstrated using one-photon absorption processes 
between quantized energy levels within a single Joseph- 
son potential well Recently, Wallraff et al. pre- 
sented experimental evidence of multiphoton absorption 
between quantized energy levels within the single poten- 
tial well formed by a large current-biased Josephson junc- 
tion ^l^J. In this Letter, we report the first observation 
of muliphoton transitions between superposition states 
of macroscopically distinct states [l5j |. which are formed 
in the double- well potential system of a superconducting 
flux qubit. The presence of an energy gap at the degen- 
eracy point of the macroscopically distinct states will be 
shown to be a prerequisite if we are to excite the qubit 
by applying a microwave even when the operating point 
is far from the degeneracy point. 

Our device is fabricated by lithographic techniques 
that deflne the structure of an inner aluminum loop form- 
ing the qubit and an outer enclosing SQUID loop for 
the read-out (see Fig. H^a)). They are spatially sepa- 
rated but magnetically coupled by the mutual inductance 
Af ~7 pH. The inner loop contains three Josephson junc- 
tions, one of which has an area /3(=0.7) times smaller 
than the nominally identical areas of the other two. The 
larger junctions have a Josephson energy of Ej = hlc/2e, 
where Ic is the critical current of the junction and e is the 
electron charge. The outer loop contains two Josephson 
junctions. By carefully designing the junction param- 
eters 0, 1^, the inner loop can be made to behave as 
an effective TSS 0. In fact, the read-out result of the 
qubit changes greatly with the qubit design ranging from 
the purely classical to the quantum regime It is de- 
scribed by the Hamiltonian H^h = ^{eoaz + A&x), where 
ax,z are the Pauli spin operators. Two eigenstates of (Jz 
are localized states | |) corresponding to the clockwise 
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FIG. 1; (a) Scanning-electron- microscope picture of a su- 
perconducting flux qubit system, (b) Signal-to-noise ratio 
A/qb//std (solid curve) and average switching current Isv 
(dashed curve) as a function of applied flux through the 
qubit loop. The two arrows indicate the signal-to-noise ra- 
tios at 4>qubit/*&o=0.5 and 1.5. (c) Spectroscopic data of the 
qubit. Resonant frequency is plotted as a function of the half 
distance between the resonant peak and dip in $qubit/'&o. 
The solid curve represents a numerical fit to the data. The 
dashed line is the energy difference between localized states 
I i) and I t). From this fit, we obtained Ej/h=S80 [GHz] and 
A/27r=0.56 [GHz]. 
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FIG. 2: Applied magnetic flux dependence of the qubit signal 
SIsw We subtracted a sinusoidal background signal from the 
averaged switching current of the dc-SQUID 7sw. (a) Without 
microwave irradiation. The data were obtained after averag- 
ing 2000 measurements, (b) With microwave irradiation. The 
data were obtained after averaging 500 measurements. The 
microwave power and frequency were 5 dBm at the generator 
and 9.1 GHz, respectively. Resonant peaks of up to two- 
photon processes and dips of up to three-photon processes 
were observed. 



persistent current of the qubit and | f ) corresponding to 
the counterclockwise current. Energy eigenstates |0) and 
1) of iJqb show energy anticrossing with energy gap HA. 
An externally applied static magnetic flux <&qubit gener- 
ates the energy imbalance heo = ^p$o(^qubit/<f o — /op) 
between the two potential wells, where $o = h/2e is 
the flux quantum, /p — Ic\/l — (l/2a)^ is the persistent 
current of the qubit, and /op is a qubit operating point 
which is a half integer. The effective energy gap in the bi- 
ased qubit is hAh — 'h\J Sq + A^. The qubit dynamics is 
controlled by a time-dependent RF field s{t) = s cos uj^xt 
provided by an on-chip RF-linc. This leads to an addi- 
tional term in the Hamiltonian i/Rpl^) = ~^s{t)<Tz- If 
A is zero, i/qb + Hrf(^) has only diagonal terms. In 
this case, the RF field cannot excite the qubit from the 
ground to the first excited state. A finite A is there- 
fore a prerequisite for a spectroscopic experiment. In 
other words, resonant peaks and dips in the qubit signals 
are direct evidence for the coherent superposition of the 
macroscopically distinct states and for the multiphoton 
transitions to occur between them. 

The qubit state was detected by measuring the switch- 
ing currents of the dc-SQUID. We defined the switching 
current as the current when the I-V characteristic ex- 
ceeded 30 ^V. The probability distribution of the switch- 
ing current was obtained by repeating the measurements 
typically 500 times. The measurements were carried out 



in a dilution refrigerator at a temperature of 30 mK. 
The bias voltage had a triangular waveform and was fed 
through large bias resistors of 1 Mil to achieve a current- 
bias measurement. The sweep rate of the current was set 
at 120 /iA/s and the frequency of the wave was 310 Hz. 
The RF line used to apply microwaves contained three at- 
tenuators in the input line: 30 dB at room temperature, 
10 dB at 4.2 K and 10 dB at the base temperature. We 
fabricated an on-chip strip line to achieve strong coupling 
between the qubit and the RF line. 

To achieve a good read-out resolution, we chose the 
operating point fop at <I>qubit/^o=l-5. Figure^b) shows 
the signal-to- noise ratio (S/N) A/qb//std and the av- 
erage switching current Isw as functions of the applied 
flux through the qubit loop, where /std denotes the stan- 
dard deviation of the SQUID switching currents over 150 
events, which may be considered to be a noise level in the 
qubit readout. The qubit signal amplitude A/qb should 
be proportional to the fiux derivative of Isw, and this 
means that we can estimate A/qb by using the signal 
amplitude AJqbo at $qubit/*i'o=l-5 (see Fig. |ll[a)). In 
experiments, the qubit signal appears when <&qubit/$o is 
a half-integer and the two arrows in Fig.QJb) indicate the 
S/N at <I>qubit/^o=0-5 and 1.5. Hence we chose the oper- 
ating point at 4'qubit/*I'o=l-5 to achieve a higher readout 
resolution. 
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Figure [2Ia) shows the qubit signal SIsw as a function 
of the external flux <i>qubit at <I>qubit/^o — 1-5, which is 
derived by subtracting the sinusoidal background from 
the SQUID switching current Jgw It shows a change in 
the thermal averaged persistent current of the qubit. We 
obtained the sample temperature of 65 mK from a nu- 
merical fit to the data by using the qubit parameters Ej 
and A, which are derived from the spectroscopy mea- 
surements (see Fig. n^c)). Figure Elb) shows the $qubit 
dependence of the qubit signal under microwave irradi- 
ation. The observed distinct resonant peaks and dips 
are attributed to situations, in which the effective energy 
separation of the two qubit states Ab matches an inte- 
ger multiple of the RF photon energy nhcucx- We have 
detected up to three resonant peaks and dips for various 
fixed RF frequencies. It should be noted that, at the 
resonant peaks and dips, the qubit is in a macroscopic 
quantum superposition of the two energy eigenstates. 

The half width at half maximum (IIWIIM„) and 
the normalized amplitude An of the dips are shown in 
Fig. |2Ia) and (b) with various microwave amplitude /rf 
for n=l, 2, 3, which is defined by /rf = exp[PR,F/20]. 
Here Prf [dBm] is the microwave power at the signal 
generator. We derived the IIWIIM„ and An from fitting 
the $qubit/*I'o dependence of <5/sw at around the reso- 
nant dips (see Fig. Elb)) by using a Lorentzian with a 
linear background. The An values are normalized by the 
fuU amplitude of the dips (3 [nA]). The HWHM„ car- 
ries important information on the dephasing and relax- 
ation processes caused by the interaction with the envi- 
ronment. To calculate the line shape of the one-photon 
absorption dip, we have used phenomenological Bloch 
equations in the first approach. To describe the strong 
driving region and the n-th dip, we replace the single- 
photon Rabi frequency in the Bloch equations by the cor- 
responding multi-photon frequency derived from theory 
of the dressed-atom approach . This substitution has 
been verified by numerical simulations ^ 19j . In a second 
approach, we have used a microscopic driven spin-boson 
model with a standard Ohmic spectral density '3,'23| to 
calculate the line-shape of the resonances. Both meth- 
ods yield similar results in excellent agreement with the 
experimental ones as seen in Fig.|3Ja). 

The Bloch equations describe the dynamics of a spin 
i in a constant field in the z-direction and a time- 
dependent field perpendicular to it. To keep the anal- 
ogy between the spin ^ and the flux qubit, we write the 
Hamiltonian in terms of the energy eigenstates and 
obtain 



H = 



h 



Ab - j^sit) 
Ab 



(1) 



The second term of the Hamiltonian leads to a non- 
adiabatic periodic variation in the Larmor frequency. 
When — Ab, the Rabi frequency is decreased. How- 
ever, we can disregard this term under the usual exper- 
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FIG. 3: Half width at half maximum (a) and normalized am- 
plitude (b) of the resonant dips as functions of the microwave 
amplitude 7rf. The microwave frequency is iL!x/2n = 3.8 
[GHz]. Solid (dashed) curves represent theoretical fits ob- 
tained by the Bloch equations combined with the dressed- 
atom approach (by the real-time path-integral method). We 
obtained rri=31(27) [ns], Tr2=84(66) [ns], r,3=430(240) [ns], 
r^=4.0(4.2) [ns], c=0.55(0.55). We have obtained similar re- 
suhs for cj^/27r = 9.1 [GHz] and 11.4 [GHz] (not shown). 



imental condition, -^s < Ab 



We have confirmed this 
by numerical simulations. If we adopt the rotating- wave 
approximation, the motion of the qubit spin (fT(i)) in the 
laboratory frame can be described by the Bloch equations 



dt 
dt 



= [-i{a{t))xB{t)],/y 
= [^{d{t))-KB{t%- 



(t)) 



((7^(<)) - CTo 



(2) 



Here, jB{t) = -^{^s cosucj, sincjcx*, Ab), cto is 
the thermal equilibrium value of {az(t)), and Tj- and 
are the relaxation and dephasing times. The steady-state 
solutions of Eqs. (HJ can be obtained in the rotating frame 
by setting '^'^"^Jl^^^ = 0. A resonant dip with a Lorentzian 
line shape appears at around Wex — Ab. The HWHMi 
and the amplitude Ai of the resonant dip follow as 



HWHMi 

Ai 




(3) 
(4) 
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where uji = is the Rabi frequency of the one photon 
absorption process. 

To describe the regime of strong driving and n-photon 
absorption processes, we apply the dressed-atom aproach 
to the flux qubit. The Hamiltonian is given by if = 
^^qb + hujcy^a^a + g(Jz{a + a^), where a and are the 
annihilation and creation operators for a field mode 
with angular frequency Wox and 5 is a coupling con- 
stant between the qubit and the field. The Hamiltonian 
can be explicitly diagonalized when A=0. The eigen- 
states are given by | T (i); iV)drosscd = exp[-(+)5(at - 
a)/?ia;cx]| T and their eigenenergies are Ejj^'' — 

NhuJex — ff^/?iWcx + {—)^heo, where |A^) are eigenstates 
of hcue^a^a. For A ^ Eq, first-order perturbation the- 
ory in the qubit Hamiltonian yields the term linear in 
A and when nhue^ ~ ?iAb, the two dressed states show 
anticrossing due to off-diagonal couphng ^ J„(a). Here, 
Jn{oL) is the n-th order Bessel function of first kind and 
a = 4:g^J {N)/?iLOex is the scaled amplitude of the driving 
field. If we prepare a localized initial state, Rabi oscilla- 
tions occur with frequency w„ = |A| J„(a) [2^ . 

We also applied a second approach starting from a mi- 
croscopic driven spin-boson model with weak coupling 
to an Ohmic bath U . In particular, we simplify the 
high-frequency approximation by numerically solving the 
pole equation (8) of Ref. [lOl- We find that the solu- 
tions for the n— photon transition are given hy 9 = 9q 
and 9 — 9n- Inserting this in the expression for Poo, 
we obtain an expression for the HWHM„ similar to that 
in Eq. (PJ), where the phenomenological dephasing and 
relaxation times T^/r follow as the inverse of the weak- 
coupling rates of the spin-boson model This result 
differs slightly from that of the Bloch equations, which 
can be recovered by setting the field-dressed level spacing 
Ao = A Jo (a) « A in the second approach. 

To analyze HWHM„ and An of the resonant dips, we 
substituted w„ for Eqs. iPJ and Q). We took Tri, T^2, 
Tr3, Tcjyn — T^, and c as fitting parameters. Here Tm 
and T^n are the relaxation and dephasing times related 
to the n-th photon absorption process. We defined a 
coupling constant c as cJrf = en. We found excellent 
agreement with the experimental data for both HWHM„ 
and An (see Fig. 13). It should be noted that the relax- 
ation and dephasing times and do not depend on 
the frequency of the corresponding multi-photon tran- 
sition for a pure structure-less Ohmic environment (see 
also Eqs. 10). This is no longer the case for a more 
complicated structured harmonic environment jl^ as it 
is present in our device. Here, the plasmon frequency of 
the dc-SQUID provides an additional energy scale of the 
environment. Nevertheless, the global physics is captured 
by our simplified model. 

In conclusion, we have reported on measurements of 
macroscopic superconducting circuits that reveal their 



quantum-mechanical behavior at low temperature. We 
have identified multiphoton transition processes in the 
qubit and found that the width of the n-photon reso- 
nance scales with the n-th Bessel function with its argu- 
ment given as the ratio of the driving-field strength to 
the frequency of the photons. Our results add another 
example to the few cases of systems that exhibit clear 
quantum-mechanical behavior at a macroscopic scale. 
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